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We report on calculations of broadening effects in QCL due to alloy scattering. The output of nu- 
merical calculations of alloy broadened Landau levels compare favorably with calculations performed 
at the self-consistent Born approximation. Results for Landau level width and optical absorption 
are presented. A disorder activated forbidden transition becomes significant in the vicinity of cross- 
ings of Landau levels which belong to different subbands. A study of the time dependent survival 
probability in the lowest Landau level of the excited subband is performed. It is shown that at 
resonance the population relaxation occurs in a subpicosecond scale. 

PACS numbers: 73.21. La, 73.21.Cd, 78.67.Hc 



I. INTRODUCTION 

A salient feature of the Quantum Cascade Lasers 
(QCL) is the modulation of their light output by a mag- 
netic field applied parallel to the growth axisi^^. This 
modulation can even be so strong that the laser turns 
off, a phenomenon refereed to as QCL blinking. The 
magnetic field Landau quantizes the electron in - plane 
motion and the laser output shows a resonant decrease 
whenever the ground Landau level of the upper state of 
the lasing transition becomes energy lined up with an ex- 
cited LL of the lower (or intermediate) state or with one 
LO phonon replica of such an excited LL. The scattering 
effects between LL's of a 2D gas are singular because of 
the unperturbed LL's are evenly spaced and macroscop- 
ically degenerate. They have been mostly studied inside 
the ground subband for degenerate electrons (the situa- 
tion of the Quantum Hall Effect, see e.g.—). The blinking 
QCL's are particularly interesting because the scatter- 
ings produce spectacular effects and involve a change in 
the subband index and of a LL index for a non degen- 
erate electron population, fn this paper, we report on 
the calculations of the inter LL, intersubband scatter- 
ing effects due to alloy fluctuations. The alloy scattering 
is known to be a very efficient scattering mechanism in 
quantum wells 7 and it is of importance to ascertain its ef- 
fect for mid infrared QCL's. Both a numerical study and 
a semi - analytical study will be performed. In section II, 
we discuss our model as well as the numerical tools we 
have used in the numerical calculations. Section III will 
be devoted to results for the density of states and the 
LL widths while in section IV we shall present results 
on the optical absorption associated with intersubband 
transitions. In section V, we discuss the time dependent 
survival probability of an electron in the ground LL of 
the upper state. 



II. MODEL 

We work in the effective mass approximation and want 
to solve the Schrodinger equation 



H*(r,t) = ih-^{r,t) (I) 
at 

with 

+ V coni (z) + KlloyW 

where we have used the transverse gauge for the vector 
potential and neglected spin effects. The z dependent 
problem admits several subbands 7 - : 

(j^ + V con{ (zf)xi(z) = E lX i(z) (3) 

For mid infrared QCL's, the subband energy distance 
is ~ O.IeV. Thus, it is a good approximation to as- 
sume that the alloy scattering is going to destroy nei- 
ther the subband structure nor at high field the LL 
structure. Thus we can consider a basis that includes 
a few subbands and LL's. Moreover, we shall replace 
in the following the elaborate layer sequence of an ac- 
tual QCL by a single quantum well with infinite barriers. 
The well thickness is L z . The numerical computation 
is done by diagonalizing the Hamiltonian in a large box 
(lOOnm x lOOnm x L z ). To model the Gai-^In^As alloy, 
we partition the large box into tiny cubes (5 A side). In 
each cube the scattering potential is a random variable 
equal to xAV with probability 1 — x and to —(I — x)AV 
with probability x where AV = 0.6eV— : 

Klloy = 5>AVf (V-i? Ga ) 

-Y,(l-x)AVf(?-R In ) (4) 

-Ri„ 

where f(Ri) is a step function which is one in the tiny 
cube centered at a lattice site Ri and zero elsewhere. 
On the scale of the Landau levels or the z dependent 
subbands \n the function / is well approximated by 
LUoS(f— Ri) where lvq is the volume of the unit cell. There 
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is no correlation between the x values of different cubes. 
The computer generates a random occupation of the tiny 
cubes. An average over N such trials is done at the end 
of the calculations. N was taken equal to 100. A numeri- 
cal diagonalization is used to extract the eigenvalues and 
eigenf unctions of the perturbed LL's. The calculations 
are done with m* = 0.05m , x = 0.53, E 2 —Ei = 143meV 
and L z = 125A In the following we concentrate on the 
vicinity of the critical field B r that lines up the (2, 0) LL 
and the (1,2) LL (we label (l,n) the n th LL related to 
the I th subband). With the material parameters we use, 
this field is B r = 30.9T. 



III. DENSITY OF STATES AND WIDTH OF 
BROADENED LANDAU LEVELS 

On the scale of the cyclotron radius or well width, the 
alloy fluctuations act like delta scatterers. With the Lan- 
dau gauge, the unperturbed LL eigenstates of the Z-th z 
dependent subband are given by: 
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where H n is the n-th Hermite function. Each of the 
(l,n) LL carries a degeneracy 2 * x % and has an unper- 
turbed energy e^ n . In our confining box this means 48 
states at 20T. In the presence of short range scatterers 
and provided one restrict the analysis to intra-n (and a 
fortiori intra-Z) transitions, one gets at the self consis- 
tent Born approximation (SCBA; see e.g£) the averaged 
Green function and self energy : 



Pi,n(e) = - 



2tt 2 HV 



(7) 
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and extends over ~Vi n ,in on each side of the unper- 
turbed (l,n) LL energy. The matrix element Vi n ,in av- 
eraged over the position of the alloy fluctuation is such 
that: 
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As it is well known Vi n ,in does not depend on n. Note 
that in our case of an infinitely deep quantum well, Vi n ^i n 
does not depend on I either since: 



dz xf(z) = for all I 



(9) 



The assumption of this model, namely that the impu- 
rity potential has only intra-n, intra-1 matrix elements is 
reasonable in strong magnetic field for the usual 2D gases. 
In the case of QCL's, where one is interested in the fate 
of electrons staying on the n — LL of an excited sub- 
band (called here E 2 for the sake of dcfiniteness), this 
approximation fails when this n = LL is intersecting 
one LL (here n = 2) that belongs to a different subband 
(called here E\ for the sake of definiteness) . 

To analyze this situation, one can restrict the Hilbert 
space to these two close-by LL's and write the wave func- 
tion of a state: 
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£+ - £i, n - S/.„(e) + iO 
Si,„(e) - vl hln (G(e)) ln (6) 



C V (ky)(pi t 2,ky + d U (ky)ip2,0,k v (l0) 



where c v {k y ) and d v {k y ) are constant coefficients and 
v labels the eigenstates. If we restrict ourselves to the 
SCBA we have now to solve: 



vf , 

v Ln,ln 
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where ( ) avcr stands for the average over the alloy fluc- 
tuations. Hence : 
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The averaged density of states (DOS) of the (I, n) 
broadened LL retains the familiar semi-elliptical shape 9 
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These coupled self consistent equations can be decou- 
pled and transformed into quartic equations for either 
£2,0 or Si,2- They are cumbersome to solve except in 
special cases; for instance when one looks exactly the 
resonance field B r where the two unperturbed LL's cross 
i.e. Ex t 2 — £2,0- If, in addition, we look at the center of 
the broadened LL's e — £1,2 = £2,0, we get : 



£2,0 (e) 



Si,a (£) 



v 2 



— 2 

2(1. 'in ^20.12 



£2,0 (e) Si, 2 (£) 
— 2 — 2 

^12,12 ^12,20 



(13) 



Si,2 (e) S 2 , (g) 



where: 
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The density of states at this particular energy s is equal 
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us with the eigenstates, eigenenergies and therefore DOS. 
We show on the upper and lower panels of fig[T] the nu- 
merically calculated DOS for non resonant fields (energy 
distances between LL's ^> LL widths) for different LL's 
of different subbands. We see that they are almost identi- 
cal, which confirms the reliability of the SCBA in the case 
of isolated LL's: not only is the DOS shape independent 
of the LL index but it is found to be also independent of 
the subband index. 

In the case of coupling and at resonance (B = B ri mid- 
dle panel in fig. [1]), we find that the LL's have acquired 
an extra width due to intersubband alloy scattering. In- 
deed, as shown on fig[U the height of the DOS peak is 
not twice that of a single peak. The numerical ratio be- 
tween the peak height of the coupled LL's (rhs peak) to 
that of a single LL (lhs peak) is ~ 1.53. This value is very 
close from the SCBA result (2-^/3/5) and again witnesses 
the reliability of that approach for short range scatterers, 
even including resonant couplings between LL's. 



IV. INTER-SUBBAND 
MAGNETO-ABSORPTION BETWEEN ALLOY 
BROADENED LANDAU LEVELS 



Pt{e) = p 2 ,o(£) + Pi,2® = 5=5* ,_ 2 2 _ 2 (15) 

'V ^20,20 + ^20,12 

If there were no coupling between the (2, 0) and (1, 2) 
LL's and at energy e for either LL at B = B r one would 
get: 

P2,o(£-£) = tx" ^ =pi )2 (£ = £) (16) 

Hence, for these particular values of the energy (e) 
and field (B r ) the ratio between the density of states 
of the alloy - coupled LL's to that of either of the LL's 
neglecting interaction is : 

— = — =2,/| ~ 1.55 (17) 

P2,0 Pi, 2 V 

Similarly, it can be shown that at B = B r , the DOS 

retains a semi - elliptical shape and extends over on each 
/ 2 2 

sides of 2y 1/20,20 + ^20,12- We note that like in the iso- 
lated LL case the coupled LL width at resonance varies 
like y/B^. When we numerically compute A, the full 
width at 1/e maximum, we find a A(B = B r = 30. 9T) ~ 
9.6meV. Note that if there were only intra-(2, 0) contri- 
butions to the width of the Landau level, one would get 
A ~ 7.5meV while if only inter (2,0) — (1,2) contribu- 
tions were retained this width would be 6.1meV. This 
demonstrates the efficiency of the alloy coupling between 
LL's belonging to different subbands and warns against 
an oversimplified addition of the scattering contributions 
to LL width (or collision frequencies). 

A numerical solution of the Schrodinger equation ei- 
ther inside a single LL or by taking into account the 
coupling between the two (2,0) and (1,2) LL's provides 



The optical absorption probability from the (1,0) LL 
states to the (2,0) LL is closely related to the density 
of states and directly measurable. In the absence of dis- 
order, the only allowed optical transition connects one 
n = LL state of E\ to the same LL state of E 2 - It 
is allowed for light propagating in the layer plane, i. c. 
in the z polarization ((X1MX2) 7^ 0). Since there is no 
non-parabolicity in the LL energies, the optical absorp- 
tion probability spectrum a(u) of an ideal sample is a 
delta function : 

a{uo) = AS (fiuj - (E 2 - £1)) (18) 

We show in figH] the absorption lineshape calculated 
according to: 

a(w)oc5^ \( i &v\z\* n =Q,v')\ 2 

x5{e v -e v/ -&«;) (19) 

where v labels the final eigenstates displayed in eqflOl 
while v' labels the initial states \& n =o.u') of the alloy- 
broadened (1,0) LL. 

In figdl we see that the absorption lineshape in the 
vicinity of the E% — E± energy resembles a broadened 
delta function. Its full width at half maximum is shown 
in fig [3] versus B where we see that increases roughly like 
\B when the two interacting levels are sufficiently energy 
separated. It is worth remarking that the DOS of each LL 
of the initial and final states has roughly the same width 
off resonance, in agreement with the SCBA predictions. 
Instead, in the vicinity of the resonance the absorption 
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width acquires an extra contribution (about 1.5meV) due 
to the inter-level coupling induced by the alloy disorder. 
The latter also activates optical transitions between the 
(1,0) and the (1,2) levels that are normally forbidden 
in ideal structures. This cyclotron resonance harmonics 
occurs in the wrong polarization (e eBl / /z) and involves a 
An = 2 forbidden change in the LL index. Fig. [3] shows 
that this disorder - induced transition gains a substantial 
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Figure 1: Density of states for different values of the magnetic 
field : B = 18 F (top), B = B r = 30.9 T (center) and 
B — 35 F. (bottom). Solid lines correspond to the DOS 
considering both the (2,0) and (1,2) LL states, dashed lines 
are associated to the DOS with only (1,0) LL states. For 
B = f8 F and B = 35 F, we show the (l,n) index of the 
states that mainly contribute to each DOS peak. Away from 
the resonance field B r (top and bottom figures), the shape of 
DOS depends neither on the LL nor on the subband. 



strength in the vicinity of the B r field. 



V. TIME DEPENDENT SURVIVAL 
PROBABILITY IN THE (2, 0) LANDAU LEVEL 

One of the key issues in QCL physics is the quantita- 
tive understanding of the escape of the carriers from the 
upper state of the lasing transition. In a simplified rate 
equation approach, before lasing, the (2,0) population 
7120 is obtained by solving: 
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Figure 2: Absorption for different values of the magnetic field 
: B = 22 F (top), B = B r = 30.9 F (center) and B = 35 F. 
(bottom) 
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Figure 3: Absorption width as a function of the magnetic 
field. Fhe solid line is a guide for the eyes. 
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The spontaneous emission power W is found propor- 
tional to: 



(21) 



where r r et r nr are the population relaxation times 
of the upper level of the transition due to radiative and 
non radiative phenomena respectively. Around B ri the 
radiative lifetime is not expected to vary much with an 
external magnetic held. On the other hand, the escape 
time due to non radiative processes involves elastic or in- 
elastic scatterings from the (2, 0) LL to other levels (here 
(1, n)). Thus, the oscillatory output power of the QCL's 
subjected to a strong field is a natural consequence of the 
modulation of the electronic DOS by the field. Implicit 
in the previous analysis is the assumption that once a 
carrier has left (2,0) it does not come back to (2,0) be- 
cause, once in (1, 2) (or whatever level) it will escape this 
level to a third one,..., much more quickly than return- 
ing to (2,0). In other words, it is assumed that the dy- 
namics of the intra-subband relaxation (1,2) to (1,1) for 
instance is always faster that the dynamics of the inter- 
subband transitions. This is often justified, in particular 
at zero field, because the density of states are continuous 
and the inter-subband matrix elements are consistently 
smaller than the intra-subband ones. In QCL's, one faces 
instead a situation where the DOS presents gaps and it 
may turn out that an inter-subband channel opens while 
the intra-subband ones are non existent In order to quan- 
titatively analyze this aspect of the alloy - assisted escape 
processes, we reformulate the problem by projecting the 
2 LL's hamiltonians as well as the alloy scattering po- 
tential on the 2 subspaces spanned by the 2 unperturbed 
LL's: 



eigenstate |Zo) very quickly drops to zero. But this does 
not necessarily imply that the output of the QCL will be 
much affected if most of the transitions amount to keep- 
ing the state inside the T>2,o (Ho + V a u oy ) V 2 ,o subspace. 
A good measure of the population relaxation time in the 
upper state of the lasing action is therefore provided by 
the calculation of the survival probability P s to find the 
system not in a given particular initial state \Iq) but in 
the whole V 2 . (H + Kiioy) Vi$ subspace: 



p s (t) =E*K«))I 2 =E|( / I e ~ ii " t/ '%)| ( 23 ) 

1 

FigfJ] shows the time dependence of P s (t) when \Iq) is 
the central level of the broadened (2, 0) LL's for B = 
18T, 30.9T and 35T. These three fields correspond to off 
resonant situations (18T and 35T) and the situation of 
exact resonance (30. 9T) respectively. In the case of non 
resonant situations (fig. 21 upper and lower panels) one 
finds fast beating between 1 (obtained at t = 0) and a 
value very close to 1: in practice there will not be any 
significant transfer between the two kinds of LL's and the 
laser will be bright if the only cause of depopulation of 
(2, 0) is the transfer assisted by alloy disorder to (1, 2). In 
fact, this behaviour is reminiscent of the time evolution 
of 1 level |0) coupled to N degenerate levels \ri). The 
survival probability in |0) is given by: 



Po(t) 
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where 5 is the detuning between |0) and the N degen- 
erate levels and fiil an effective coupling between |0) and 
these levels which is given by: 



m = .f£(0\V\n)(n\V\0) 



(25) 



H = Ho + V^oy (22) 

= V 2 ,0 (H + Kdloy) V 2 ,0 + Vl.2 (Hq + V^lloy) Vl,2 

where V 2 ,o and "Pi, 2 are the projectors on the unper- 
turbed (2,0) LL and (1,2) LL respectively. The two 
first terms represent the intra-LL broadening due to alloy 
scattering while the last two terms induce transitions be- 
tween the (broadened) LL's We assume that at t = the 
system has been placed in one of the eigenstates \Iq) of 
^2,0 (Ho + Kiiioy) P 2 ,o- As time flows, the system leaves 
this eigenstate. While the intra-(2, 0) disorder is included 
in the basis of the eigenstates |Z) of V 2 ,o (Ho + V^noy) V2,o> 
second (and higher) order couplings between these eigen- 
states exist due to the virtual jump(s) to the (1,2) sub- 
space. Thus, the survival probability in any particular 



As applied to the QCL's, this simple modeling can 
work only if the detuning is large compared to the width 
of the LL's. In fact, it correctly predicts that the product 
of the squared oscillation frequency by their amplitude 
should be proportional to B. In addition, it shows that 
the period of the oscillations decreases with increasing de- 
tuning between the two LL's, in agreement with the nu- 
merical calculations. However, the oscillations between 
1 and 5 2 / (5 2 + 47i 2 fi 2 ) in eqUHare everlasting while the 
numerical calculations exhibit a fast decay of the P s (t) 
envelope to a constant 1 - (2h 2 fl 2 ) / (S 2 + 4h 2 fl 2 ). This 
decay is clearly associated with the different Bohr fre- 
quencies of the problem: both the initial LL's have a 
finite width, thereby leading to many Bohr frequencies 
and thus to the replacement of the oscillatory sine func- 
tion by an exponential decay. At resonance (figHl middle 
panel), the survival probability shows a very fast decay 
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(subpicosecond) to a value 1/2. At large time, the sur- 
vival probability in the initial subspace is just the fraction 
of states that belong to this subspace. In our calculations 
there 74 states in each subset, thus the 1/2 limiting value. 

We have checked that the previous conclusions remain 
valid when the initial state is not the central state. Actu- 
ally, whatever the initial state, the calculated time con- 
stants (pseudo-period, characteristic transfer time at res- 
onance) are independent of the initial state within a few 
%. This proves that our findings regarding a given ini- 
tial state can be extended to the population of the ini- 
tial Landau level irrespective of its precise nature (e.g. 
Boltzmann-like) . In actual samples, there are of course 
other mechanisms that empty (1, 2) but do not refill the 
(2, 0) LL and finally, the survival probability in the ini- 
tial subspace will of course goes to zero. The key point 
is that the system leaves the initial subspace in a very 
short time. Whether, in actual samples, this escape is 
irreversible or partly irreversible (which would occur if 
the escape from Ei, n = 2 is too slow) remains to be 
elucidated. In the case of a quasi reversible escape, like 
fast elastic scatterings between LL's, we believe the pop- 



ulation relaxation is better described by: 

Tr V P+Ah 2 W) ^ ' 

where S = £2,0 — -^1,2 and HQ is an effective coupling 
constant between the (2, 0) and (1, 2) LL's and more gen- 
erally between (2,0) and (l,n). A minimum is clearly 
reached at resonance, but this minimum is non zero. To 
suppress the QCL emission, it is necessary to include 
other losses. A preliminary version of this work was pre- 
sented at the ICPS 28 Vienna (2006). 
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Figure 4: Survival probability P 3 of states for different values 
of the magnetic field : B = 18 T (top), B = B r — 30.9 T 
(center) and B = 35 T. (bottom). The initial state is taken 
at the central level of the broadened (2, 0) LL's 



